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ABSTRACT 



The problem of finding the optimal waterjet propulsion 
pump for a hydrofoil of given configuration and displace- 
ment is examined from the point of view of least weight. 
Included in the propulsion pump study are the gear weight 
for matching the waterjet pump to the specified prime 
mover (s) and fuel weight for a specified range and speed. 

From the current waterjet pump designs, two competing 
design approaches were selected: (1) Multi-parallel double 

suction impeller centrifugal pump with up to ten impellers 
and (2). Axial pump with an inducer impeller followed by 
axial stages when required by the pump head and net positive 
suction head. Because of the hydrofoil's drag to speed 
characteristics with its take off hump, the pump design 
problem requires an efficient reliable pump that meets 
the cavitation conditions at take off. Both designs 
satisfy these requirements. 

From the required flow rate and system losses each 
pump design is solved for size and weight with the choice 
of optimum pump being the lightest design for that particu- 
lar set of overall system parameters. Each pump design is 
directly related to the cavitation condition at the take 
off point which is the design point. The value of Thoma's 
cavitation criterion is closely related to the pump selec- 
tion. As the value decreases the selection of pump shifts 
from the centrifugal to the axial design. The actual 
selection for a given design point is a trade off between 
the pump, gear, and fuel weights where the centrifugal 
design has the heavier fuel and lighter gear weight. 

Thesis Supervisor :A . Douglas Carmichael 

Title: Professor of Ocean Engineering 
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Chapter 1 INTRODUCTION 

Waterjet propulsion is not a new ship propulsor and 
has been used or tested on a variety of ship types ranging 
from the three thousand ton naval destroyer to the small 
pleasure planing craft. Besides pleasure craft, waterjets 
have been most successfully used by river patrol craft in 
South Vietnam. 

The large ocean going hydrofoil has brought renewed 
interest to waterjet propulsion because of the unique 
propulsion problems presented by the hydrofoil. Propeller 
propulsion with the power transmission through the strut 
has posed many problems to the designer and is currently 
considered to be unreliable. The waterjet has been 
demonstrated as a feasible and reliable system on the 
hydrofoil Tucumcari (PGH-2) . Low overall efficiency and 
high plant weight are still problems to the system. 

The current contract by the U.S. Navy for a waterjet 
hydrofoil in the two hundred ton range has caused increased 
interest and development work by industry on waterjet 
systems. Also, considerable work has been done on water- 
jets for surface effect ships which also encounter the 
same transmission problems and high craft speeds. 

To the pump designer the waterjet system requires a 
departure from the normal practice of pump design and has 
made adapting off-the-shelf pumps to hydrofoils a difficult 
process. It has also created entirely new designs or 
caused adaptation of radical pump types from other fields. 
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Because of its low specific weight and small size, 
the aircraft derivative gas turbine has become the 
principle prime mover for the hydrofoil in its foil 
born mode. Even when British engines are included, the 
selection of gas turbines with respect to a shaft horse- 
power is very limited in the ranges from five to ten 
thousand and fifteen to twenty thousand shaft horsepower. 

The problem of matching a given prime mover to the 
optimum pump type and number of each is a complex and 
not well understood problem. The problem has been 
studied to define the design parameters and how they 
relate to give the optimum design based on least system 
weight. The weight of the propulsion package will include: 
pump, gear, and engine dry weight; ducting and piping; 
entrained water above the water line; and fuel weight for 
a' given hydrofoil range at constant craft speed. The pro- 
blem is bound by the hydrofoil configuration, size, and 
speed. The configuration will be similar to the PGH-2 with 
the two after struts containing the inlet nacelles with the 
prime mover (s) and pump(s) located internally in the craft. 
The maximum number of engines and pumps will be two per 
strut. Only subcavitating foil speeds will be examined, 
therefore, limiting the cruise speed to fifty knots. A range 
in displacement from fifty to one thousand tons will be 
considered . 

The pump design problem will be to establish the 
lightest pump, gear, and fuel weight combination for a 
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given set of water jet propulsion parameters. All of the 
system components will be from tested designs. The pump 
design will be studied in detail while the gear and prime 
mover will only be studied for their weight growth and 
relationship to the overall problem. These designs will 
be incorporated in a computer design program. This program 
will be a subroutine in an optimization program for the 
overall waterjet propulsion system. 
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Chapter 2 WATERJET PROPULSION 

The waterjet propulsion system may be considered as 
four sub-systems: inlet nacelle, ducting, pumping machinery, 

and nozzle. In order to understand the selection of pumping 
machinery, its interaction with the three other areas is 
required . 

Using the mass average form, the thrust delivered 
by the system: 

T = pQ(Vj cosa - V Q ) (2.1 ) 

where cosa is the depression angle of the nozzle. In the 
equilibrium state the thrust equals the drag for a given 
craft speed. 

The velocity ratio, Vj/V Q , is the primary independent 
variable for system evaluation. Put into the momentum 
equation and neglecting the depression angle gives: 

T = pQV 0 (Vj/V Q -1) (2.2 ) 

Figure (1) taken from Levy's report-'- is a typical study of 
propulsive efficiency versus the jet velocity ratio for 
values of system loss coefficient where the losses are 
approximated by: 

AH l = K d Vj 2 /2g (2.3 ) 

The studies show a flattening out of the efficiency curve 
to the right of the maximum efficiency for a given loss 
coefficient. This fact leads to a jet velocity higher than 
that at maximum efficiency because of the large decrease in 
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entrained water weight for a small decrease in system 
efficiency. This curve illustrates why systems designed 
on maximum efficiency tend to be heavy. 

The shape of the hydrofoil drag versus speed curve, 
see figure (2), illustrates one of the basic design 
problems. In its displacement or hull born mode, the 
hydrofoil drag increases with speed to a maximum value 
just prior to flying. This point will be referred to 
as take off and typically represents the maximum value 
of drag. Just after take off the drag drops off and then 
increases with craft speed. For the typical subcavitating 
foil system, the take off speed is about one half the foil 
born cruise speed, ie . , for a 45 knot cruise speed, take 
off speed is 22 knots. 

To allow for acceleration, the available thrust at take 
off is assumed to be 25 percent greater than the equilibrium 
thrust at the take off point. 

The main system parameters are selected outside the pump 
design problem. They set the required flow rate and system 
losses. From the parameters and losses the required pump 
head can be calculated: 

H = Vj 2 /2g - V Q 2 /2g + AH L + h^ (2.4 ) 

Also the net positive suction head, which is the head above 
vapor pressure, at the pump inlet is found for use in 
cavitation considerations for the pump. It is given by: 

H sv = Vo 2 / 2 ^ " AH L + H a " h £ -h v (2.5 ) 
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It is these values of flow rate, pump head, and net 
positive suction head along the drag curve that drive 
the pump design. 
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Chapter 3. WATERJET PUMP DESIGN 
Section 3.1 Basic Theory 



The pumping machinery suitable for water jet propulsion 
are turbomachines since other types are too heavy. Turbo- 
machinery covers a large range of pumps which can be 
classified by the flow path through the pump with respect 
to the rotating axis and grouped into three types: 
centrifugal, mixed flow, and axial. 

The principle by which the pump transmits energy 
to the fluid is the conservation of angular momentum. 

In the case of steady flow, the change in net flux of 
angular momentum through the rotating impeller is equal 
to the torque, t, on the impeller. Referring to the 
velocity diagram. Figure (3) , : 



= C pr 2 V 2 Cosa 2 dQ-/ p 1 V 1 cosajdQ (3.1 ) 

' 2 J 1 x 



For uniform flow across the impeller this reduces to 
the form: 



t = pQ(r 2 V 2 cosa 2 - ^Vjcosaj) (3.2 ) 

The power required by the pump is: 

P = Taj = p QgH (3.3 ) 

Combining equations (3.2) and (3.3) gives Euler's equation 
for the theoretical head, Hth, developed by an impeller 
without hydraulic losses: 



Hth = 1/g (wr 2 V 2 cosa 2 - ur jV jcosaj) 



19 



Hth = 1/g (U 2 V 02 - U i V Q j ) 



(3.4 ) 



Euler's equation further reduces to the form: 

Hth = 1/g (U 2 V 02 ) (3.5 ) 

when there is no prewhirl to the flow at the impeller inlet. 
This is a good approximation for most pump impellers with 
no inlet guide vanes and a uniform flow profile. 

From equation (3.5), the theoretical head to capacity 
curve can be produced which gives an idea of the pump 
characteristics with respect to discharge angle from the 
impeller with the assumption that fluid angle, B 2 , leaving 
the rotor remains constant: 

V9 2 = U 2 - V m2 cot $ 2 

Hth = U 2 2 /g (1-V m2 /U 2 cot e 2 ) (3.6 ) 

The flow rate is proportional to V 2 m . This gives a linear 

relation of head to flow rate. 

The definitions of velocity diagram angles varies 

considerably in the published literature. Figure (3) is 

2 

based on definitions used by Stepanoff and will be carried 
throughout . 

For example, a centrifugal impeller with back sloping 
blades, where 6 2 is less than 90°, and constant pump speed, 
produces a theoretical head proportional to Q: 

Hth « Cj - C 2 Q cot 6 2 ' (3.7 ) 
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This gives a negative slope straight line. The power to 
capacity can also be shown: 



P « Hth Q 




(3.8) 



Equation (3.8) gives a convex curve with some maximum power 
point. This fact could explain why back sloping impellers 
are often used. 

Euler's equation for head also leads to definitions 
of pump efficiency. Efficiencies are divided into three 
groups: Hydraulic, volumetric, and mechanical. Only the 

hydraulic efficiency effects the pump head and is defined 
as the ratio of pump head to theoretical head: 



The hydraulic losses are caused by skin friction along the 
direction of flow and eddies. The prediction of hydraulic 
losses is usually not accurate for a new design. 

Volumetric efficiency is a measure of the leakage 
losses or losses in capacity. 

Mechanical efficiency is the measure of losses in power 
that are not related to the head. They include disk fric- 
tion and bearing losses. The total pump efficiency is the 
product of the three groups: 



n h = H/Hth 



(3.9 ) 




(3.10) 



n 



gpQH (English units) 



(3.11) 



550 SHP 
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Section 3. 2 Principles of Similitude 

The current state of preliminary pump design is based 
on several established similarity conditions for geometri- 
cally similar pump designs. These relate two of our design 
input pumps; head, H, and capacity, Q, with pump size and 
speed. 

Relating the fluid dynamic forces by the velocity 
triangles at the same geometric point on an impeller gives 
the flow coefficient: 

<f>=V/uL=V/Ul 

Y m 'A m 'B (3.12) 

2 

where is proportional to Q/D and U is proportional 

to ND gives the similarity conditions: 

3 

K = Q/ND = constant (3.13) 

q 

which relates capacity to size and speed. 

From Euler's equation the forces on the pump are pro- 
portional to head: 

Hth = U 2 V 0 2 /g (3.14) 

where V0 2 aU 2 . Therefore, the head is proportional to 
tangential velocity squared. Put in non-dimensional form 
gives the head coefficient: 

<P = H /u 2 /g (3.15) 

2 2 2 

substituting N D for U gives the dynamic similarity that 

relates head to size and speed: 

K, = H/ 2 ^ 2 = constant (3.16) 

h N D 
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Since power is proportional to the product of head and 
flow rate, a similarity condition can be established from 
the above two : 

K = PA. 3 ,. 5 = constant (3.17) 

p N D 

The basic similarity conditions are only good if the 
hydraulic losses of the impellers being compared are a 
constant fraction of the pump head, i.e., no Reynolds 
number effects. In practice, the assumption is good 
especially for the range in pump size that the hydrofoil 
will be utilizing. 

The flow and head coefficients are important parameters 
in pump design. This fact will be treated in detail in the 
following section. 



Section 3.3 Specific Speed 

By combining the two similarity conditions to get a 
condition relating head, capacity, and pump speed: 



3/4 _ 



K 2/K h = nq'/h 



3/4 



(3.18) 



This parameter is specific speed, N g , and is used in 



three different sets of units: 



Ns (gpm) = N(rpm) Q (gpm) ^/H ( f t ) 3//4 
Ns(cfs) = N(rpm) Q (cfs) 2 /H(ft) 3//4 
n g = N(rps) Q (cf s) ^/ (gH) 3//4 

Due to usage, the non-dimensional form, n , is not as common 

s 

as the dimensional forms. In order to handle comparisons, 
the following conversions are given: 



(3.19) 

(3.20) 

(3.21) 
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Ns(cfs) = Ns(gpm)/21.2 



ns = Ns (gpm) /17150 

ns = Ns (cf s) /811 . 3 

Section 3.4 Pump characteristics 

Pump characteristics are given in the form of a head 
versus capacity curve and an efficiency versus capacity 
curve. Both curves are for a given pump speed. 

Since the similarity conditions, K , and, , hold 
along a constant specific speed or efficiency line, the 
pump performance map can be constructed from one speed 
line; figure (4) is given as an example. A plot of Kh 
versus Kq, figure (5), is another form of the same infor- 
mation and is sometimes more useful. 

The specific speed of a pump is usually at the point 
of best efficiency, BEP. Once the design is established 
from the operating requirements at BEP, the pump character- 
istics at off design can be found from the performance map. 

Assuming the normalized head curve with respect to BEP 
is parabolic: 

2 

H/H d = A (Q/Q ) + B Q/Q d + C (3.22) 

and converting to a Kh versus Kg curve: 

2 

K = A'K + B'K + C' (3.23) 

h q q 

the original curve becomes: ^ 24) 

2 2 
K/H d = A (Q/Q d ) + B Q/Q d N/N d + C(N/N d ) 
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In the case of the hydrofoil, knowing the design point, 
the head, and the flow rate for the off design, the off 
design pump speed can be found by solving the above 
quadratic equation. 

Another aspect of the pump performance curves is their 

general shape with respect to impeller type and specific 

speed. The hydraulic losses produce a curve rather than the 

linear relationship given by Euler's equation for head. 

Figures (6) and (7) are examples of changes in performance 

curves for a range of specific speeds for centrifugal pump 
2 

impellers . 

To explain the change in slope of the head curve at 
the best efficiency point, examine the design curves in 
figure (10) . As the specific speed increases the flow 
coefficient also increases. When the capacity, ie. , Vm , 
is increased for a constant pump speed, the percentage 
change in the tangential velocity component is much smaller 
for small values of $2 than large values of <j> 2 . Therefore, 
for high values of specific speed, the head term, which is 
proportional to the tangential velocity component, drops off 
sharply with an increase in flow rate. 

Section 3.5 Impeller Types 

Turbomachinery can be classified by three basic impel- 
ler types based on the flow path through the impeller with 
respect to the rotating axis. Also impellers are further 
classified by specific speed. Although each type could 
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theoretically cover a broad range of specific speed, usage 

has dictated certain ranges for each impeller type to 

achieve the best efficiency. These ranges are only general 

3 

and there is definitely a large overlap in each type . 

In the centrifugal impeller, the fluid exits radially 
with respect to the rotating axis by changing the flow 
direction 90 degrees through the impeller. The specific 
speed range is from 24 to 189. This indicates that the 
centrifugal pump is used for high head and lower flow rate 
applications . 

The mixed flow impeller has an exit flow path somewhere 
between the radial and axial direction. Its range is from 
189 to 472. 

As the name suggest, the axial pump flow path is 
parallel with the pump axis. Its specific speed range is 
472 to 708. The axial pump is used for higher flow rate 
and lower head conditions where additional stages are 
added to cover higher head requirements. Much of the axial 
gas turbine compressor design work has influenced the axial 
pump design. 

All three of the impeller types are represented in 
waterjet propulsion. For small craft application, Jacuzzi 
builds a pump with one mixed flow impeller and Buehler uses 
a two stage axial design. The current hydrofoil, PGH-2, 
uses a centrifugal pump with two double suction impellers 
operating in parallel on one shaft. By splitting the flow, 
the centrifugal pump can handle the large flow rates. 
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Further discussion of the pump types with their design 
parameters and performance characteristics is covered in the 
design section. 

Section 3.6 Cavitation 

In classical fluid dynamics, cavitation occurs when 
the local static pressure of the fluid is equal to the vapor 
pressure: 

2 2 (2 25) 

P + JspV = Pv + h pV 1 ' 

To the pump designer, cavitation imposes severe 
limitations. The vapor bubbles forming along the pump 
impeller can cause impeller damage by collapsing and 
impinging on the surface eroding the metal. This process 
can also cause vibration and noise in the pump. A fall off 
in pump head is caused by a more developed vapor cavity. 

The third input to the design problem, net positive 
suction head, Hsv, is the most important for cavitation 
considerations since it gives the total head above vapor 
pressure at the impeller inlet: 

2 

Hsv = h-hv + V m ! /2g (3.26) 

The point at which cavitation occurs for given types 
of impellers is determined from experimentation. The point 
where there is a three percent drop in head is most common 
and will be used here. Incipient cavitation would have 
already occured at this condition. 

Prandtl's cavitation number in terms of relative inlet 
velocity at the impeller shroud is: 
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P-Pv 



(3.27) 



a = 2 

P ^PW S 1 

Substituting in the equation for Hsv gives: 

W 2 V 2 

Hsv = - 2 ±- + K (3.28) 

2g 2g 

where is taken to be 1.1 to account for impeller eye 
conditions. This equation relates the velocity inlet 
triangle, <f> j , to Hsv for a given cavitation number. If 
is the critical value, then the minimum value of Hsv is 
determined. 

There is another cavitation constant used in pump 

2 

design that was first introduced by Thoma . Thoma said 
that there is an experimental ratio between the inlet velo- 
city head at capacity cut off and the pump head. This 
was first used for low Ns pumps. Thoma ' s constant is 
defined : 

a T = Hsv/H (3.29) 

and is constant for a given Ns for conditions approaching 
cavitation while the similarity conditions still hold. 

There is also a Reynolds number effect but some uncertainty 
exists how this effect changes in c T « A recent study 4 of 
large centrifugal impellers shows o T increasing with impel- 
ler size which is contrary to what theory suggests. 

A third parameter is found by applying the similarity 
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conditions to Hsv that were applied to pump head. This 
is a suction specific speed and relates Hsv to N, and Q: 

S = NQ^ Hsv" 3/4 (3.30) 

S is expressed in the same units as Ns with the same 
conversions holding. 

Suction specific speed is a useful parameter for it 
relates the cavitation limitations to specific speed: 

Ns = Sa T 3/4 (3.31) 

For example, if a limiting value of S is used and Thoma ' s 
criterion is determined from the pump requirements, then 
the pump specific speed is set. 

When suction specific speed is expressed in terms of 
cavitation number, inlet flow coefficient, and inlet 
geometry, as developed in appendix A, some insight is 
given into the influence of cavitation on the pump design: 

60 { w(i„ r 2))Js 

IT ** 

s = 

<f> x { l/2g (1 + o + )} j/4 (3.32) 

P 2 
♦ l 

Figure (9) is a plot of the above equation and shows how 
there is an optimum value of tf> i , ie., inlet angle, for a 
given value of cavitation number. Also, it shows that a 
minimum value of hub to tip ratio gives a maximum value 
of S. 
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5 

Wislicenus states that for commercial pump operation, 

a value of a for cavitation free operation is 0.4. However, 
P 

the Hydraulic Institute uses S equal to 8500, which corres- 
ponds to a 0 ^ of 0.2 as the limit. Specially designed pump§ 
which will be covered in detail, have reached values of S 
equal to 34,000. As figure (9) shows, these pumps have low 

values of <j> i and x and operate in the cavitation region. 

P 

Cavitation is the primary limit to hydrofoil water jet 
pump design. Furthermore, because of the increased losses 
and lower ram head, the take off point is the critical point 
for cavitation and, therefore, drives the pump design. 
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Section 3.7 Pump Design Problem 



The curve of hydrofoil drag to ship's speed establishes 
a two point design problem for the pump designer. One point 
is at takeoff and the second at cruise. The takeoff point 
requires the pump to maintain head and avoid cavitation 
damage and the cruise point requires good pump efficiency. 
There is a choice of designing at either point and checking 
the other for cavitation and efficiency. Since the inlet 
conditions are critical for establishing good performance 
against cavitation, the takeoff point is chosen as the 
design point. 

Current approaches to pump design for water jet pro- 
pulsion can be put in three catagories: 

(1) The first design is the centrifugal double 
suction impellers in parallel on one shaft. 

This design is demonstrated by the (PGH-2) 

pump which was built by Byron- Jackson . 

4 

Lockheed has also recommended this design. 

The design approach is to reduce the flow rate 
per impeller. This raises the pump speed 
and decreases the impeller size. This design 
will be covered in detail. 

(2) The second design approach is to have a two 
speed multi-stage pump. The first stage 
operates at a lower RPM which allows the 
pump to meet the inlet cavitation require- 
ment. The remaining stages have a higher 
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speed and produce the majority of pump head. 

This design is being developed by both Aerojet 
and Garrett. The disadvantages are the special 
gear arrangement and double shaft. This arran- 
gement is feasible but requires development to 
insure reliability. This design will not be 
covered since it is similar in size and weight 
to the third design. 

(3) The third approach is the single speed multi- 
stage pump with the first stage being an 
inducer which operates at a suction specific 
speed approaching 30,000. This design has 
been built and tested by Pratt and Whitney 
and by Rocketdyne. The inducer has proven 
both efficient and reliable, after early 
problems in the inducer design. 

The pump design section of the optimization program will 
calculate the design parameters, size, and weight for both 
the centrifugal, multi-parallel, double suction impeller 
design and the axial inducer design. The lightest weight 
pump including entrained water, gear, and fuel weight is 
sub-optimized within the overall program. 

Each pump design will be established at the design 
point, ie. , take off point. From the relationships of 
Thoma's cavitation criterion, specific speed, suction 
specific speed, and design coefficients to the required 
flow rate and head at take off, the pump size and weight 
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are established for each design. Any other point along the 
drag curve will be considered an off design point and the 
pump speed and efficiency are calculated from the pump 
characteristic head and efficiency curves. The following 
design sections are based on current pump designs incorpo- 
rating realistic values of suction specific speed at the 
takeoff point. 

Section 3.8 Centrifugal Pump Design 

For hydrofoil application, the centrifugal pump has 
been selected by several pump designers and manufacturers. 
They use multi-parallel, double suction impellers on one 
shaft and in one pump casing. The PGH-2, which has one 
pump with two double suction impellers, has demonstrated 
that the design is both reliable and efficient. 

The multi-parallel design has several advantages over 
the single stage pump design when both are operating near 
the same suction specific speed. Given the operating con- 
ditions, they both have the same Ns for one impeller: 

Ns = Sa T 3//4 (3.33) 



For the double suction design, the specific speed will be 
for one impeller, ie. , one half of a double impeller. With, 
n, being the number of double suction impellers, as n in- 



creases and impeller size decreases: 

3/4 



N = 



NsH 



Qh 



(2n) 



(3.34) 



This reduces gear weight and produces a long narrow pump 
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design when six to ten double suction impellers are used. 

The number of double suction impellers depends on output 
speed from the prime mover and available space and machinery 
configuration in the hydrofoil. For design purposes, a 
maximum of ten impellers is used. 

The program objective in the pump design is to esta- 
blish the pump size and weight. The pump dry weight is a 
function of the pump geometry and its size. The geometry 
reflects both the specific speed and the type of discharge 
system that is used in the design. The design is based on 
the pump model (ref. 4) that was developed specifically for 
waterjet propulsion. Since the pump configuration and dis- 
charge system (Fig. 10) is established the weight is only 
a function of size and specific speed. 

By letting the impeller exit diameter, D 2 s, represents 
the size, the dry weight is: 

V D 2s = F (N ) (3.35) 

s 

The value of x is somewhere between two and three depending 
on the growth and casing thickness. From the stress analy- 
sis in reference (4) the minimum casing thickness is much 
less than the thickness used for economical casting techni- 
ques. Therefore, it is assumed that for the range in pump 
head experience in hydrofoil design, the casing thickness 

is a constant. This fact gives x = 2. This realtionship 

1 

is further supported by the Peerless pumps plotted in 
figure(H). Even though their weight is heavier than 
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the water jet pump designs, they demonstrate that equation 
(3.35) for x - 2 is representative. 

Knowing the head coefficient of the pump establishes 
the exit diameter since pump speed is known: 

D 2s = Is- ( 3 . 36 ) 

However, the head coefficient is not constant over the range 
of specific speed but varies slightly with . By substitu- 
ting into the definition of specific speed, a relationship 
can be found between the head coefficient, flow coefficient, 
and pump geometry term at impeller exit. From appendix B: 

n = (Ka b/D ) H i> 2 h /^ z 3/A 

o 2 ^ o 

(3 37) 

where Ka is the blockage factor. 

From current pump design practice the impeller discharge 

angle is a constant in the range of 22 degrees to 23 degrees 

for good pump performance. By setting B 2 constant, the 

specific speed can be expressed as a function of 2 an( 3 

b/D where: 

4 s 

4 1 2 = n h (1 - <P 2 cot B 2 ) (3.38) 

assuming no inlet prewhirl and a value of hydraulic efficien- 
cy. The plot of equation (3.37) for values of B 2 anc ^ 
specific speed is often used as a design chart for centri- 

2 f 6 

fugal pump design . However, before this plot can be 

completed something must be known about the impeller width 

ratio, b/D , with respect to N . As illustrated in figure 
^ s s 

( 12 ), the impeller exit width and blade curvature change 



35 



considerably as N g is increased. Both the curves from the 
2 7 

literature ' and values of width ratio from current designs 
show the same slope of b/D 2s with respect to specific speed 
(Fig. 12). By establishing the width ratio the design chart 
(Fig. 13) can be completed and the head coefficient for a 
given specific speed can be found giving the impeller exit 
diameter . 

From figure (11) , the pump weight for a single section 
is established from the three current designs available: 



W n = C N D 2 (3.39) 

D w s 2s 

where C =3.32 for a single discharge volute, 
w 2 

Following current practices for light weight pump 

construction, the casing is cast from aluminum with a 

titanium impeller and steel pump shaft. From the Lockheed 
4. 



report a typical weight breakdown for an eight double 
suction impeller pump with an N = 57.2(1210): 



Pump case 


2260 


lbs. 


Impellers 


832 




Shaft 


627 




Bearings 


431 




Misc . 


350 




Total 


4500 


lbs. 



For good pump efficiency the impeller inlet angle was 
taken equal to the exit angle, ie., 8 i= 3 2 * The requirement 
for higher suction specific speed causes a smaller value 
of inlet angle to be used. In order to study the possibil- 
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ity of increasing the suction specific speed for the 
double suction impeller design, Lockheed tested two 
impellers with different configurations. Using the normal 
configuration for inlet ducting and a five bladed impeller 
with an inlet blade angle of 13 degrees and an exit angle 
of 22 degrees, a suction specific speed of 13,800 (gpm) 
with a Hsv of 31 feet was obtained. This relationship of 
low flow coefficient to high S was established in section 
(3.6). For current centrifugal design an inlet of 13 to 15 
degrees is considered a minimum. Therefore, for this de- 
sign, Bi was set at 14 degrees, <(>2 = 0.25. By knowing the 
blockage and hub to tip ratio at the inlet, the flow coef- 
ficient can be used to find the inlet diameter: 



240 Q/n 1//3 

Fair 2 <(> 1 N(l-r hl 2 ) 



(3.40) 



For the multi-parallel design the minimum hub diameter 

is established by the pump shaft requirements to handle the 

input torque. A value of r^ = 0.5 will be used which 

4 

corresponds to designs established by Lockheed . A detailed 
design would require an increase in r^ for increases in the 
number of double suction impellers. As the number of im- 
pellers increases, Q per impeller decreases, thereby 
decreasing D^ s , but the power transmitted remains constant, 
holding minimum constant. 

The centrifugal pump length is scaled from the impeller 
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inlet diameter, D^ g . The pump length is based on the 
section length for one impeller (Fig. 10) plus the shaft 
and bearing supports at the pump ends. The section length 
is governed by the inlet system. Letting the average inlet 
velocity at the pump base be equal to 75 percent of the 
impeller inlet velocity, V gives the impeller section 
inlet area: 

A = 

0.75V, 

lm 



A = 



Q/N 60 

0 . 75uN D <f> , 
ls Y 1 



(3.41) 



By letting the rectangular inlet be: 



A = L L 
x y 



(3.42) 



where : 



by letting: 



L - pump section length 



Ly - pump section width 



L = C t L 
y Lx 



where : 



'L = 0.866 



then: 



L = /A/C. 

X ' 3 



The total pump length is then: 



L = nL + L 
p xx 



(3.43) 
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The pump efficiency is then found by scaling for 
Reynolds number from a standard value of = 0.88 
when D„ =2.33 feet. The efficiency scaling law is 

4 S 

9 

based on a method used by Byron Jackson : 



^h 





D 



2s 



0.165 



(3.44) 



The centrifugal pump off design point characteristics 
are found by assuming parabolic head and efficiency curves. 
The head curve for the centrifugal design changes with 
changes in specific speed. The off design pump speed is 
found by using equation ( 3.24 ) where the coefficients 

are a function of specific speed. This relationship is 
plotted in figure (14) . 
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Section 3.9 Inducer Design 

The problem of pumping the cryogenic liquids for 
rocket engines created the need for an impeller that is 
generally referred to as an inducer. The rocket fuel pump 
was required to match a high rpm input, usually a vapor 
turbine, to handle the high flow rates and to be as small 
and light weight as possible. In order to reduce size, 
the pump speed must be increased. This requires an increase 
in inlet suction specific speed for a given pump head and 
net positive suction head at the pump inlet. To increase 
S a special inducer impeller was put in front of a centri- 
fugal impeller raising the Hsv for the centrifugal impeller. 

The same requirements for the rocket pump are also 
present for the hydrofoil water jet pump. Therefore, the 
inducer design has been adapted to the waterjet pump appli- 
cation . 

The design goal is to operate at the highest possible 
suction specific speed, ie., high rpm, without reducing 
pump efficiency or reliability. The inlet geometry and, 
therefore, pump size is governed by the equation (3.32) 
developed in appendix A. As stated in the cavitation 
section, low values of cavitation number, hub to tip ratio, 
and inlet flow coefficient are required to achieve high 
values of suction specific speed at the pump inlet. 

Figure (9) shows that there is an optimum inlet blade 
angle, ie., <f> } , for a required suction specific speed. 
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Taking the derivative of equation (3.32) with respect to 

<f> j gives the optimum value of 4> ^ for a given value of 

a , see appendix A. Achieved values of a are determined 
P P 

experimentally and depend on the amount of impeller loading 
along the blade. The ratio of maximum inlet relative 
velocity to actual relative velocity at the tip can be put 
as a function of cavitation number: 



a 

P 



Wmax 

( ) 

W 



2 



“1 



(3.45) 



For low values of a , approaches l, ie., no loading 



at the impeller inlet. For S = 30,000, a very low value of 

a is required. If r, = 0.3 and a = 0.03, then the opti- 
p hi p 

mum value of <f> j is 0.166 from equation (3.32) and a theo- 
retical suction specific speed of 26,000 should be achieved. 

From the inlet analysis, the inducer size and inlet 
geometry are fixed. The outlet is fixed by the second 
stage design. Typically, the next stage has a lower 
suction specific speed impeller and, therefore, its inlet 
flow coefficient is much higher. If the second stage has 
the same tip diameter, then the merdial velocity must be 
much larger. This causes the hub to tip ratio to increase 
by a factor of two or larger through the inducer. Some 
designs also decrease the impeller tip diameter through 
the inducer. This improves flow to the second stage and 
decreases its tip speed. This is carried out when there 
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is a large change in flow coefficient. 

The blade shape of the impeller is very critical to 
its performance. Several design methods are in the 
literature, references (10) to (16), but the best designs 

have been achieved through experimental and developmental 

10 11 
work ' 

The major problem in adapting the inducer type impeller 
for hydrofoils is in preventing cavitation damage while 
not sacrificing cavitation performance. Since a rocket 
is only a short duration vehicle, maintaining pump head 
was the primary concern and some impeller erosion was 
acceptable . 

There are several current designs by industry of 

successful uses of the inducer impeller. Sponsered by 

the Joint Surface Effect Ship Project Office, both Pratt 
11 10 

and Whitney and Rocketdyne have developed and tested 
pumps that have good cavitation performance and resistance 
to cavitation damage. Since both of these inducers have 
approximately the same minimum o and specific speed, they 
will be used as the basis of the first stage of the axial 
pump design. The impeller parameters used are: 



N 



s 



r 

r 



hi 

h2 



4i 2 



149.5 



0.3 



0.7 

0.4 
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(j) 2 ” 0.11 



a Tmin = 0.055 

Section 3.10 Axial Pump Design 

The axial pump consists of the inducer impeller 
discussed in section 3.9 followed by axial stages when 
needed (Fig. 15 and Fig. 16) . The inducer can be used in 
front of any of the three impeller types. The axial stage 
was selected because for many cases only the inducer im- 
peller is required to handle the pump head and a,p. This 
fact is due to the speed restriction imposed by the sub- 
cavitating foils. Only by going to very high jet velocity 
ratios does the pump head exceed 700 feet. The axial 
stage design is the lightest weight due to the inlet 
matching and discharge nozzle associated with the design. 
The nozzle section can be as short as 0.65 times the inlet 
diameter. These two factors give advantages over the 
centrifugal design where the discharge system is a large 
part of the pump weight and the mixed flow pump where a 
change in tip diameter and flow direction require a longer 
and more complex problem of matching the stages. 

The axial stage design consists of a stator row and 
rotor. Patterned after the axial compressor design, the 
design has a higher hub to tip ratio, a higher flow coeffi 
cient, and a lower head coefficient than the inducer stage 
A series of detailed designs and tests, references (17) 
through (19) , give three designs with hub-tip ratios of 
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0.4, 0.7, and 0.8 with increasing blade loading as the 
hub-tip ratio increases. 

The axial stage must be matched to the inducer. From 
the inducer exit hub to tip ratio and the ratio of inlet 
to exit tip diameters, the axial stage inlet flow coeffi- 
cient can be calculated: 



< ^ >1 axial C ^ 1 inducer 
2 2 

£ _{Ka D lc , (l~r^ ) linducer 

{Ka D (l-r ) }axial 

is hi 



(3.46) 



For Ka and D to be constant for both designs if * =0.11 

is * v 1 i 

then <|> = 0.196. If the tip diameter is decreased through 

the inducer as in the Rocketdyne design, the r^ is increa- 
sed for the inducer raising the flow coefficient. From 
the axial designs, the parameter diffusion factor is used 
as a measure of rotor blade loading: 



DF = 1 



W , R V - R V 
2S + S2 fl 2 SI 01 

Wi c/sW CR 



+R ) 
S2 s l 



(3.47) 



for R = R , DF reduces to: 

S2 SI 



DF 



W 2 + A V 0 
W ls 2 (c/s) Wj 



(3.48) 



As shown in fugure (17) the diffusion factor gives the 
relationship between impeller loading and flow and head 
coefficients. Because of the high head inducer design, 
a DF of around 0.5 is used giving ^ 2 = 0.2. Using esta- 
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blished values of solidity, c/s = 1.01 for the r^ = 0.7; 
the preliminary design is established from the flow and 
head coefficients. 

Figure (19) is a plot of head and flow coefficients 

versus specific speed for inducer designs available in the 
10 11 12 13 

literature ' ' ' ' . Since the higher head inducers, 

which are the Rocketdyne and Pratt and Whitney designs, 
are in a narrow range of specific speed; this specific 
speed value is used for the inducer design. By using 
one inducer the design approach to the axial pump is 
modified by letting determine the suction specific 
speed. 

Depending on Thoma ' s cavitation criterion, the take 
off suction specific speed can vary up to a maximum value. 
When the maximum S is exceeded, an additional axial stage 
is added to the pump design. 

Using flow rate and pump head, pump speed is calcula- 
ted. Inlet diameter is calculated from N and flow coeffi- 
cient . 

To determine the number of stages both o T and head 
coefficient must be used. When either the minimum value 
of or maximum value of S are exceeded, the number of 
stages increases. 

Pump length and dry weight are a function of the 
number of stages and inlet diameter: 

L =C t D (3.49) 

P L si 
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X 



(3.50) 



W 



Dt 



CD 

W 



is 



where x is less than three and greater than two. Pratt & 

l 

Whitney sets X = 2.2 and other sources set X at 2.5 for 
more conventional designs. Since a series of pump sizes 
is not available for this design, a value of X = 2.3 will 
be used. Scaling from the Rocketdyne pump weight, the 
coefficients are: 



No. Stages 



C 

w 



C 



L 



1 347.1 1.79 

2 393.5 1.71 

3 439.9 2.03 



The off design operating points are calculated by 
applying the similarity conditions to the characteristic 
head curve. Figure (4) is the map for the inducer design 
and was generated from equation (3.24). From the off 
design rpm, the ratio of off design to design flow coef- 
ficients is calculated. If the flow coefficient ratio is 
much less than 1.0, stall is possible; and if it is greater 
than 1.0, cavitation is more likely or a higher Hsv is 
required. Since off design points are at a higher craft 
speed, ie., increased ram head, the off design cavitation 
problem does not occur. From test cases run, there has 
not been a significant change in cj> j from the take off to 
cruise operating point. 

Pump efficiency based on measured results is set at 
91% for a 3.66 foot impeller diameter. The efficiency is 
scaled for larger and smaller pump sizes using equation 
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(3.44). The off design efficiency is found by assuming 
the efficiency curve with respect to flow rate to be 
parabolic . 
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Chapter 4 REDUCTION GEAR DESIGN 

The lightest weight gear design is the gas turbine's 
integral gear box, and therefore, when available it should 
be used. When a reduction gear is required to match the 
pump speed to the prime mover output shaft speed, a gear 
design based more on aircraft than marine practice is 
designed . 

The purpose of the gear design is to give the weight 
growth associated with changes in pump speed and different 
engine-pump combinations. 

Three different gear designs are incorporated in the 
program to handle the possible range in gear ratio and the 
various engine to pump combinations. When the number of 
engines equals the number of pumps, a planetary gear is 
incorporated because it has the lightest weight. For an 
unequal number of pumps and engines a combining gear is 
required to match either two engines to one pump or one 
engine to two pumps. This gear is based on a single reduc- 
tion offset gear with idler design for gear reduction ratios 
less than twelve and a double reduction double branch design 
for higher gear ratios. 

The gear weight determination is based on a method 
20 

developed by Willis . The method uses the assumptions 

that the gear weight is proportional to the solid rotor 
2 

volume, Fd , and the rotor volume is a function of input 
torque, gear ratio, and K factor where K factor is uniform 
for all gears: 
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SHP and output speed, N^, the gear weight is calculated by: 



w = o 35 126,000 SHP ( £Fd j 

G ■ ~ 1 - — - - — - — - - 



K 



Np 



(4.5) 



For the planetary offset gear designs, the K factor 
equals 500. This requires a hardness of = 60 which 
implies a hardened and ground gear. In order to have a 
light weight and relatively small gear, the high K factor 
is essential. To demonstrate the feasibility of a large 
planetary gear, Curtiss Wright built a 40,000 SHP gear for 
the U.S. Navy. 

The offset combining gear design is based on a single 
input and output shaft design times a constant to account 
for the added pinion and idler gear in the case of the two 
engine one pump combination or the added bull gear for the 
two pump one engine combination. For the two pinion gear 
design the weight would be less if the total SHP 
were used because the power is divided decreasing the pinion 
and idler size. However, the value of SHP used in the 
program is always SHP per engine. Therefore, the weight 
increases by a factor of 1.3. For the two bull gear design 
the weight is 1.7 times that of the single shaft design. 

The double reduction gear design is included in order 
to cover entire range of gear ratio. Current marine locked 
train double reduction gears are too heavy for hydrofoil 
application. Therefore, a light weight gear design of 
higher K factors requires development. Because the K 
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3 



(4.1 ) 



2 

d F 



x/2k 




For each gear type an expression of the overall gear ratio 
M , and the first reduction or idler gear ratio when 

applicable is developed and set equal to the weight factor 

2 2 



EFd = EFd 
x/ 2 k C 



(4.2 ) 



An example of the design method is the single offset gear 

with idler where the weight factor is equal to: 

2 2 

EFd _ , 1 + nu + mu + Mo + Mo 

— i g g 

C ra g m g 

(4.3 ) 



To find the lightest gear design the value of m^ is deter- 
mined by taking the derivative of the right hand side of 
equation (4.3) with respect to m^ and set equal to zero. 

For the offset design this reduces to: 

3 2 2 

2mg + rr g = Mo + 1 (4.4 ) 

2 

Similar expressions are developed for all gear types and 
are summarized in table (2). 

From the weight factor, gear weight is found by 
determining the desired K factor which is a measure of 
compressive tooth stress and corresponds to the surface 
hardness. An allocation factor, AF, is also used to adapt 
the general design method to particular applications. 

For hydrofoils, AF = 0.35. Therefore, for a given engine 
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factor used in the design method is an average value and 
also to bias the problem against the larger gear, a lower 
K factor value of 300 is used. This value would still 
require a harden gear for the first reduction pinion and 
gear . 

For the planetary and offset gears, the gear efficien- 
cy is set at 0.98. For the double reduction gear =0.95. 

Figure (19) illustrates the influence of gear ratio 
on gear weight for the four designs. The gear weight 
divided by input torque is plotted versus the gear ratio. 
From the figure it is clear that the planetary gear is 
by far the lightest design. 
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Chapter 5 PRIME MOVERS 

The aircraft derivative gas turbine is the predominant 
prime mover for hydrofoil craft. Table (1) is a list of 
British and American engines that have been marinized and 
are, therefore, suitable for hydrofoil propulsion. The 
list is not all-inclusive but gives the range of power 
available. For example, the Rolls Royce Olympus is not 
included but the FT4 and LM2500 engines cover that parti- 
cular power range. 

A confusing aspect of comparing gas turbines is the 
different power ratings. For this problem, only the max- 
imum intermittent power and the normal power ratings are 
used. Another commonly used term is maximum continuous 
power which lies between the two listed ratings. The inlet 
ambient temperature plus inlet and exhaust losses also 
effect the power ratings. The three common ambient tempera- 
tures used by manufactures are 60° F, 80° F, and 100° F. 

The table reflects the current United States Navy practice 
using 100° F inlet temperature as the standard. 

Some engines are available with an integral reduction 
gear. The Proteus engine is listed for the 1500 rpm and 
1000 rpm versions of output shaft speed. 

The matching of gas turbine to craft displacement is 
critical for an optimum propulsion plant weight to displace- 
ment ratio. Because of the existing gaps in available en- 
gine sizes, ideal matching of engine power to required power 
is not possible for all displacements. Instead of determi- 
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ning craft size by mission and required payload, the avail- 
ability of a gas turbine or combination of gas turbines 
should be a primary consideration in the decision of craft 
size. 

The engine size, efficiency and output rpm all have 
an impact on the overall design. However, the single 
largest weight component that the engine influences is fuel 
weight. Depending on the required hydrofoil range, the 
fuel weight can be a significant percentage of the overall 
propulsion plant weight. The fuel weight is calculated for 
a constant cruise speed assuming that the lift to drag 
ratio is constant and the system losses are only a function 
of the jet velocity head term. 

The measure of engine efficiency is given in terms of 
specific fuel consumption, SFC. A characteristic of the 
gas turbine's performance is that its minimum SFC is at 
maximum power. The SFC then drops off as the power de- 
creases. For the fuel calculation, the relationship of SFC 
to SHP is approximated by: 

SFC = C FS /(SIIP) n (5.1 ) 

where n = h for SHP from normal power to 70 percent of 
normal and n = 3/4 for the range from 70 percent downward. 
Figure (20) shows this approximation with respect to an 
actual gas turbine curve. 

The fuel consumption for constant SHP for a time At is 

W_ = SFC SHP At 

r 
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(1-n) 



At 



(5.2 ) 



VL 



C 



FS 



However, for the case of constantly decreasing SHP as the 
fuel is burned, the fuel weight solution is not straight 
forward. The required SHP at any time is a function of 
displacement and system losses if D/L and V q are held 
constant. As developed in Appendix E: 



SHP 



pg Q H 
550 n n 

p g 



(5.3 ) 



SHP = 



yA.V . 

3 3 

550 n n 

P g 



V . v 



(h^ + (1+K^) — - 



(5.4 ) 



2g 2g 



where : 



V j = “ V o + r V o + A (D/L) 1 

2 1 4 pA. 



h 



(5.5 ) 



For solving by computer, a discrete step approach was used. 
The total fuel was then the sum of the intervals: 



W 



F 



N 

I 

i=l 



W. 



Fi 



(5.6 ) 



For a range of 750 n. miles, the difference in fuel for 
N = 20 or N = 40 was less than one half of a percent, 
therefore, 20 intervals is used for the calculations. 

As would be expected if the shaft horsepower per prime 
mover is less than 70% of normal power at full load, at the 
end of 750 n. miles the required SHP will only be 50% of 
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normal power for a craft speed of 45 knots. As figure (20) 
demonstrates, the SFC in the 70% to 50% range is much in- 
creased. These relationships will be further discussed in 
the chapter on integrated system results. 

The engine output shaft speed is important for gear 
weight. The lightest gear is the integral box and is 
treated as part of the engine weight. 

The last area of engine influence in design is the 
number of prime movers. The design method only treats 
the cases of 1,2, or 4 engines. When the number of en- 
gines increases, provided a sacrifice in SFC is not requi- 
red, the design weight tends to improve mainly due to the 
decrease in pump weight and water weight. However, multi- 
ple engines also offer advantages in flexibility and 
reliability. For example, if one engine out of four could 
be shut down during cruise, then the SHP per engine of the 
other three would increase and total engine hours would 
decrease . 

The detailed design problems of engine and pump 
installation including auxiliary systems have not been 
treated. Since the particular engine choice is outside 
of the problem, its auxiliary weights will not change. 

A detailed study of the best engine and pump arrangement 
given the number and type would prove worthwhile for future 
study. 
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Chapter 6 INTEGRATED SYSTEMS RESULTS 

The waterjet system is most strongly influenced by 
the choice of gas turbine for a particular hydrofoil dis- 
placement. As already discussed in chapter 5, the propul- 
sion system weight to displacement ratio varies considera- 
bly with displacement primarily due to the availability of 
gas turbines over a broad range in power. 

In order to analyze the design program results and 
its choice of optimum pump and design parameters, it must 
be remembered that the craft displacement and choice of 
gas turbine are inputs to the program and, therefore, are 
not optimized. Most of the analysis of the design program 
was carried out on a hydrofoil of 400 ton displacement 
because of the availability of published figures of estima- 
ted losses and weight data for that particular size hydro- 
foil. Computations were also made for 80 tons, 200 tons, 
and 750 tons of displacement. 

The program selects the lightest pump, gear, and fuel 
weight combination for a given operating condition. 

Although the size and performance of the nacelle and struct 
diffusion are important in defining the operating condi- 
tions, the jet velocity ratio gives the clearest indication 
of what type pump will be selected and what weight is 
driving the optimization program. As figure 1 illustrates 
there is a value of V^/V^ for maximum propulsive efficiency 
and values to the right of this point decrease the propul- 
sive efficiency thereby increasing the required power. 
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As stated earlier, the optimum design point will be 
for higher values of Vi/V^ than the maximum efficiency 
point. The results support this fact except when the 
problem is bounded. From table (3), it is seen that the 
400 ton hydrofoil went to a high optimum value of V j/ V Q 
(2.83) for the case of the FT4C-2 gas turbine where there 
was no power limitation. For the other 400 ton cases using 
one engine, the upper horsepower limit at take off limited 
the V^/V q to some value less than 2.83. If two engines 
were used splitting the horsepower so that only 50% of 
normal power would be required at cruise, the SFC to SUP 
relationship will drive the V^/V^ value very high as in the 




3.80. 

The pump selection can also be related to the jet 
velocity ratio when no bounds are present. The losses 
increase with V./V ; therefore, the value of Thoma 1 s 

J O 

cavitation criterion, o^, decreases. As a general rule for 
high values of a T (greater than . l ) the centrifugal pump is 

selected. For low values of o^,(less than 0.05) the axial 
design is selected. In the region between these values of 
o^, the pump size, gear weight, and fuel weight or problem 
bounds can enter into the selection process. 

Since the primary objective of multi-parallel design 
is to increase pump speed, the rpm for the centrifugal pump 
is higher than the axial design and, therefore, the gear 
weight can be the predominant factor in selecting the 
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centrifugal pump design. Pump speed can also be a bound 
for small hydrofoil displacements when the axial pump, 
due to its fixed specific speed, cannot achieve a gear 
ratio of one or greater and is eliminated. 

The fuel weight component for the same V^/V^ is a 
function of pump efficiency. The axial design is generally 
three percent more efficient, and therefore, the axial 
pump requires less fuel than the centrifugal pump. Depen- 
ding on the range, the fuel weight can also be a determining 
factor in the pump selection. 

The comparison of the pump weight component is a 
function of a and pump size. Because the centrifugal 
pump specific speed is proportional to o^, the difference 
in pump weights decreases as a T decreases until a cross 
over point is reached. 

The relationship in pump diameter also affects the 
difference in pump weight. For the same operating condi- 
tions, if the number of pumps is increased, the axial pump 
becomes lighter when compared to the centrifugal design. 
However, this decrease in weight of the axial pump can be 
offset by increases in gear weight especially if the 
increase in pumps cause a shift from a planetary to a 
combining reduction gear. 

The problem bounds are quite often the determining 
factors in setting the Vi /V q and therefore in selecting 
the pump type. The primary bound is available engine horse- 
power. Depending on the relationship between maximum 
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intermittent power and normal power, either maximum 
intermittent power can be a limit at take off or normal 
power can be a limit at cruise. It is even possible for 
the difference in pump efficiency to be a determining 
factor in selecting the axial pump for high values of a T 
because the centrifugal design is rejected due to required 
cruise power. 

The cross over region between the pump designs can 
be changed by adjusting the values of maximum suction 
specific speed and the maximum number of centrifugal 
impellers. For example, if the maximum suction specific 
speed for the centrifugal design is increased to 10,000; 
the specific speed would increase for the same values 
of a^, and the cross over point from the centrifugal to 
axial pump will lower. 

As table 3 illustrates there is a large variation 
in the design parameters at the optimum point for different 
displacement and engine selections. There are several 
explanations for this variation. First is the maximum 
available horsepower, which limits the jet velocity ratio. 
By increasing the jet velocity, the required power is 
increased. If not enough power is available, the design 
will reach a maximum value of jet velocity ratio before 
an optimum point is reached. 

The second case is the installed engine power greatly 
exceeds the full load cruise power. As in the example of 
the 200 ton hydrofoil with an LM2500 gas turbine installed, 
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the optimum point is driven to a very high value of jet 
velocity ratio typically greater than three. This fact 
is explained by referring to the gas turbine performance 
characteristics, figure (20). When the required power 
at full load is less than 70% of the normal engine power, 
the SFC to SHP relationship is such that lower fuel weights 
can be achieved at higher values of SHP, ie., Vj/V Q . 

Also, the high jet velocity ratio gives low entrained 
water weight. If the range were reduced or set to zero, 
a very different optimum point would be selected by the 
program. 

In summary, the choice of the lightest pump design 
is directly related to the optimum point parameters that 
the program selects. This selection is a very complex 
trade off between the system components. Only the cases 
that can be directly attributed to the pump or prime mover 
have been dealt with here. However, from data points 
already run, if no power limitation exists, the value of 
jet velocity ratio at the optimum design point is greater 
than the maximum propulsive efficiency point and usually 
higher than 2.0. 

The selection of pump type is primarily dependent on 
the system losses, ie., Thoma's cavitation criterion, a . 
Where for very low values of adless than 0.0 5) only the 
inducer design is selected due to the poor cavitation con- 
ditions which cause the centrifugal design to go to very 
small values of specific speed. However, as the cavitation 
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conditions improve, the centrifugal design becomes more 
competative and is usually selected for values of o 
greater than 1.0. 

In the middle range of o T near the cross over point, 
the pump selection is more arbitrary. Other factors of the 
design problem such as arrangement of machinery and weight 
distribution should then enter into the pump selection 
process . 
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Chapter 7. CONCLUSIONS 

From the study of v/aterjet pumps for hydrofoil propul- 
sion several conclusions to the design approach can be made. 
The first conclusion is that the take off point cavitation 
condition establishes the pump design. If only a conven- 
tional impeller design, with respect to suction specific 
speed, is used, the specific speed for one impeller is 
established by the required purnp head and net positive 
suction head at the pump inlet. From sample cases run for 
S = 424.5(9000), the specific speed range was from 
49 to 118. This specific speed range specifies a contri- 
fugal type impeller especially at the lower end of the 
range where the optimum point is more likely to be esta- 
blished. To design for a single impeller pump at the low 
required value of specific speed would give a large diameter 
slow speed pump v/hich would require a large reduction gear. 

In order to get a light weight pump design, the three 
design approaches stated in section 3.1 have been tried by 
industry. All three design approaches are an attempt to 
increase the pump speed and, therefore, decrease the pump 
size and weight. 

The selection of the optimum pump design from the two 
selected design approaches is very closely related to 
Thoma 1 s cavitation criterion, o , as discussed in chapter 7. 
Changes or improvements to either or both of the pump de- 
signs would only shift the cross over point and should not 
alter the overall outcome of the propulsion pump problem. 
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Chapter 8 • RECOMMENDATIONS 

The water jet pump problem has only been approached 
from the point of view of least weight for a given hydro- 
foil displacement. It would be beneficial to study water- 
jet propulsion if the problem were approached from the 
view point of lightest propulsion system weight for a 
required payload and the optimum ' propulsion system and 
hydrofoil displacement for a specified gas turbine. 

Since the program only solves the preliminary design 
for two design approaches, the updating of the current 
program as more developed work becomes available and the 
inclusion of the two speed pump design approach if develop- 
ment proves the design feasible should be accomplished. In 
the jet velocity ratio range where the two designs are 
nearly equal in weight, a detailed pump design study 
concentration on the engine and pump arrangement, location 
of jet nozzle, and details of the auxiliary propulsion 
systems is needed. 

As stated, the pump selection is limited by the choice 
of gas turbine. Only in the case of proper matching of 
engine power to displacement is a true optimum design 
achieved. Any further work should approach the problem 
from the view point of gas turbine selection and further 
study the problem when a good matching is not possible. 
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The centrifugal multi-parallel impeller design always 
requires less gear reduction than the inducer design. 

When a combining gear is required, this fact produces 
the centrifugal pump as the optimum design over a larger 
range of a T » 
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Integral R. Gear output 



Table 2. Reduction Gear Weight Equations 



Gear Type 


Gear Ratio Equation mg 


Simple Offset 


2m 3 + m 2 = 1 

g g 


Offset with 
Idler 


2m 3 +m 2 =M 2 +l 
g g o 


Double-Reduction 

Double-Branch 


2m 2 M 2 + 1 

2m 3 + 9 - ° 

g M +1 . M +1 

(0)2(0) 

1 M 0 M o 


Planetary 


t 3 2 0.4 (Mo — 1 ) 2 + 1 

2m J + m z — 

s s b 


Gear Type 


2 

Total Weight Fraction ZFd/c = 


Simple Offset 


1 + l/m + ro + m 2 

g g g 


Offset with 
Idler 


1 + l/m + m + m 2 + M 2 /m +M 2 

g g g o g o 


Double-Reduction 

Double-Branch 


h + l/2m + 2m + m 2 + m 2 /M + 

g g g g o 

M 2 /2m + M /2 

o g o 


Planetary 


1/b + 1/bm + m + m 2 + 
s s s 

0.4 (M - l) 2 0.4 (M - l) 2 

° + ° 
b m b 

s 
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able 3 (continued) 
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429 910 27 .030 149.5 1265 4.75 4.79 8 681 725 .917 .919 



Table 3 (continued) 
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Figure 1. Propulsive Efficiency versus Jet Velocity Ratio 
for Values of Loss Coefficient 
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Firgure 3 





Impeller Flow Velocity Diagram 
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% Flov/ Rate, Q/Q D 
Figure 4 . Inducer Stage Performance Map 
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Figure 5. K as a Function of K for the Axial Pump Design 



Head, H/H 




Figure 6. Centrifugal Pump Flow Rate vs Head for Changes in 
Specific Speed 
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Figure 7. Centrifugal Pump Efficiency Curves for Various 
Values of Specific Speed 
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Figure 8. Pump Impeller Profiles for Range in Specific Speed 
from the Hydraulic Institute 
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Specific Speed 
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W/D- /Impeller Section, lb/ft 



X - Peerless 
CD - PGM- 8 4 
A - PGH-2 

O - 550 ton Hydrofoil 




10 2 4 6 8 10 2 2 

Specific Speed, N g (cfs) 

Figure 11. Centrifugal Pump Section Dry Weight to Exit 



Diameter Ratio as a Function of Specific Speed 
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Figure 13. Centrifugal Impeller Design Chart for b/D =0.001N -0.025 




Coefficients, A,B,C 



Figure 14. Centrifugal Pump Head Curve Coefficients as a 
Function of Non-Dimensional Specific Speed 
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Figure 15. Inducer Stage Axial Pump Design 
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Figure 16. Inducer Plus One Axial Stage, Axial Pump Design 




Figure 17. Theoretical Head Coefficient vs Flow 
Coefficient for Value of Diffusion Factor 
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Coefficient as a Function on Non-Dimensional Specific 
Speed for Various Designs 
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Figure 20. Percent SHP vs Percent SFC with Respect to Normal 



Power Approximation for Gas Turbine 
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Figure 19. Gear Weight to Input Torque Ratio as a Function 



of Reduction Gear Ratio for FAC = 1.0 
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Appendix 

( 1 .) 



^ • Development Of Cavitation Parameters 

Prandlt's Cavitation Number Is: 

- P ~ P, 



° l = 



v 






when related to the impeller inlet 



a = ^1 ^v 

p — 






Is 



(A.l) 



a = h - hv 
P 



W ls/ 2g 



relating 
Hsv = h 



a 

P 



to Hsv 

2 



hv + 



Vml 

2g 



H, 2 V ,2 

Hsv = a ls + K1 ml 
P 



2g 



2g 



(A. 2) 



where is an allowance for induced prewhirl 

at the inlet and is typically (1.1). 

Substituting for Vml to give Hsv in terms of <}> x 

— u 



TT ,U, 2 + V . 2, , „ V , 2 

Hsv = a ( 1 ml ) + K, ml 

p 1 



2 g 



2g 



Hsv = V ml^ {a ( 1 + l/^ 2 ) + 1 

TT p 



(A. 3) 



for Ki=1.0 
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(2.) Relating S to G p and <|> j : 



% 

N Q 2 



Hsv 



J/4 



N = 60 IL/itD, 

1 Is 



Q = Kau/4 Dj (1-r^) 



S(cfs) = 60 U 1 /irD ls { Kair/4 (1-r^) 

(V ml /2g tod+l/*! 2 ) + 1 >)3/4 



S(cfs) = Ka^60/2tr^ {l-r h ^}^ 

4> ! ( l/2g ( l+o p +o p / 4> ! 2 ) } 3/4 



(A. 4) 



k 



S 1 = S/(l-r h 2) 2 



C = Ka^ 3 0 /n ^ =16.5 (for Ka = 0.95) 



, 2 3 /4 

S 1 = C/<j) 1 {l/2g(l+o p +o p /(J) 1 )} ' 



(A. 5) 



(3.) For maximum suction specific speed given a 
constant a 



I*! 4 / 3 (l+cr p +o p /<f> i )} 3/44 = minimum 

d/d<^ 1 {(f» 1 4/3 (l+a p +o p /(}) 1 2 ) 3/4 = 0 



4/3*! (1+a ) ~2/3a p 



= 0 



*1 = {Js(a p /1 +a p )}‘ 



(A. 6) 
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1 



for a - 0.03, K.. = 

P 1 

<J>! = 0.120 

for a = 0.03, K n = 1.1 
P 1 

4> ! = 0.116 
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Appendix B. Development Of Specific Speed As A Function 
Of Head and Flow Coefficients. 



(1.) Non-dimensional specific speed is defined: 

n s = n Q i ' 2 /(gH) 3//4 (B.l) 

n s = N (cfs)/811.3 

(2.) Exit area of impeller equals: 

A = Ka it b D 2s (B.2) 

where Ka is the blockage factor 
(3.) The pump speed equals: 

n = U 2 /ttD 2s (B. 3) 

(4.) Flow coefficient is defined at exit as: 

*2 = V m2 /U 2 (B.4) 

(5.) Head coefficient is defined as: 

ip 2 = gH/U 2 2 (B . 5) 

(6.) Substituting equation B.2 and B.3 into B.l 
gives : 

n s = U/»D 2s (V m2 Ka « b D 2 e )*V< 9 H) 3/4 

n s = l/» 1 (Kab/D 2s ) 1 U 2 v m2 ! / (gH) ' 

n s = l/it* 5 (Kab/D 2s )* s (V ml /U 2 ) V (gH/U 2 2 ) 3/4 

(7.) Substituting for definitions of B.4 and B.5: 

o / a 

n = 1 /tt 2 (Kab/D„ ) 2 <p 2 V’f'2 (B.6) 

S ^ o 
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( 5 .) 



The gas turbine SFC to SHP relationship is 
approximated by: 

SFC = C Fs /SHP n (E. 6) 

where : 



n = 1/4 

for 100% to 70% of normal power 
n = 3/4 

for less than 70% of normal power. 
(6.) Substituting E.6 for E.5 gives: 



W„ = C 



Fs 



SHP (1 n) t 



(E.7) 



(7.) Using a discrete step approach where the time, 
t, is equal to the range over speed: 

At = R 



V N 
o 

N 

W p = Z 
F i=l 



W 



fi 



(E.8) 



where : 



W,, . = C„ SHP (1 n) At 
f x Fs 



(E.9) 



for V. = 
3 



-V q /2 + { V q 2 /2+ ( A-Wf (i _i) ) D / L } ^ 

{ : > 



PA. 



(E. 10) 
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Appendix C„ Sample Calculation For The Centrifugal 
Pump Design. 



(1.) Input to test routine: 
Displacement: 

Cruise Speed: 

Range: 

Configuration : 

Engine : 

Jet Velocity Ratio: 



100 tons 
45 Kts 
750 n. mi. 

2 pumps, 2 engines 
Tyne RM-1A 
2.0 



Input to 


the pump solution 


from the test 




routine : 


Cruise 


Takeoff 




Q 


123.42 


126.89 


cfs 


H 


304.66 


393.80 


ft 


Hsv 


91.02 


22.56 


ft 



(3.) Pump design solution at takeoff: 
° T = Hs/H 
a T = 22.56/393.80 
o T = 0.057 
N s (cfs) = S(cfs) 



for S 



max 



424.5 

N (cfs) 
s 



49.72 (for single suction) 
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for the number of double suction impellers n=10 



Q/IMP = 126.89/2-10-2 
Q/IMP = 3.172 cfs 

N = N s H 3//4 / Q/IMP* 2 



N = 2467 rpm 

for values of <(>2 & i> 2 enter the design chart 

for = 0.90, the chart is based on b/D 2s = 

0.001 N - 0.025 for 3 2 = 22°: 
s 



<j) 2 = 0 . 16 



\p 2 = 0.552 



D 



2s 



D 



2s 



60_ ,gH. h 

* N ^ p 

60 . g 398.8, 2 

ir24 65 { 0.552 ’ 



for : 



D n = 1.172 ft 
2s 

n 240 Q/IMP , % 

Is 1 1 r 2 Ka<() 1 (l-r hl ^')'N J 

r hl = °* 5 ' 3l = 14 °' Ka = 0,95 

<() ^ = tan 14 — 0.25 

D, = 0.558 ft 
Is 



n 



P 



l-{ ( 



2.333, 

1.172 



0.165 

( 1 - 0 . 88 ) } 



n = 0.866 

P 
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SHP/ENG = 


p g Q h 

550 ri (Wo. Eng) 
1? 


SIIP/ENG : 


= 3354.931 Shp 



If takeoff SHP is less than maximum inter- 
mittent power then calculation is continued: 



Ain 


= Q 




Vin (no. pump) 


Vin 


=0.75 Vml 


Vin 


= (o.75) (0.25) ttD 1s N/60 


Ain/Pump 


= 109 Q/ND^ s (No. Pump) 


Ain/Pump 


= (109) (126.89) 

(2467) (.558) (2) 


Ain/Pump 


= 5.03 ft 2 


L 

P 


=( Ain ) (2+n) 

n 0.866 


L 

P 


= 6.85 ft 


W D 


= (0 . 725n + 0.275) C D 2 N 

w 2 s s 


W D 


= (7.525) (4.66) (1.17) 2 (49.72) 


W D 


= 2395.1 lbs. 



W = 0.55 W n 
w D 

W = 1317.3 lbs. 
w 
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(4.) Off design calculations at cruise: 



h/h d = a < QA ) d > 2 + B Q/q d n/n d + c ( n / n d ) 2 

n/n d = Ic Q/Q d ± (§c<Q/Q d > 2 - |<0 /q d ) 2 + gS^ 



A = -1.5 

C = 2 . 85n -0.17 
s 

B = 1 + 1.5 -C 
N/N d = 0.904 
N = 2225 rpm 

n p = n pD (A(Q/Q d ) 2 + B Q/Q d + C) 
A = -0.75, B = 1.51, C = 0.24 



SHP/ENG 

SHP/ENG 



0.865 

P g Q H 

550 n p (No. Eng) 

2530.39 shp 
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Sample Calculations for Axial Pump Design. 



Appendix D. 



(1.) Input to test routine: 
Displacement : 

Cruise Speed: 

Range : 

Configuration : 

Engine : 

Jet Velocity Ratio: 



100 tons 
45 Kts 
750 n. mi. 

2 pumps , 2 engines 
Tyne R1I-1A 
2.0 



( 2 . ) Input to 
routine : 

Q 

H 

Hsv 



the pump solution 

Cruise 

123.42 

304.66 

91.02 



from the test 

Takeoff 

126.89 cfs 
393.80 ft 
22.56 ft 



(3.) Pump design solution at takeoff: 



N (cfs) 
s 

Q/pump 

N 

N 

N 

M 

o 



H /H 
s / 

0.057 

149.5 

63.447 cfs 
N H 3/4 (Q/p) 15 

o 

(149.5) (398.8) 3/4 /(63.45) 3s 

1659.17 rpm 

N /N 
P 
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M = 1.874 
o 

for <f> j = 0.11, = 0.3 



D 1 

Is 


2 40 Q/pump 
~ ( } 
Tr 2 Ka(J) 1 N(l-r h 2) 


D ! 

Is 


= 2.131 ft 


Ain 


= °- 785 D ls (1 - r hl > 


Ain 


2 

= 3.275 ft 


L P 


= C T D 
Xj Is 


L 

P 


= (1.79) (2.131) 


L 

P 


= 3.81 ft 


n p 


r r 0.165 

= ( 2 ; 13 ) (1-0.915) 


n p 


= 0.907 


SHP/ENG 


_ P g Q H 
550 r) z 
p 


SHP/ENG 


= 3271.7 


SHP/ENG 


< Maximum Intermittent Power 


W 

D 


2 3 

= C D. 
w Is 

2.3 


W D 


= 347.1 (2.131) 


W D 


= 1477.98 lbs. 


W 

w 


= 0.523 Ain L pg 
P 


W 


= 414.28 lbs. 



w 
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(4.) Off design calculations at cruise for 



A = -3.0, B = 4.8, C = -0.8 
N = 1528.16 rpm 

ip = > (A (Q/ V 2 + B Q/ % + c 

for A = -1.7, B = 3.42, C = -0.72 

n =0.907 
P 

+ l'*D = ^ 1 *° 60 
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Appendix E. 



Fuel Weight Calculation Method 



(1.) 


SHP 


= pgQH/550 n p n g 


(E.l) 


(2.) 


assuming 


D/L to be constant: 






(D/L) A 


= p Q (Vj-V Q ) , ( A-lbs) 






A (D/L) 


= p A.V. (V.-V ) 
D 1 1 o' 






as fuel 
D/L, Aj, 


is burned V . is reduced with 
1 

V remain constant, 
o 


A for 




V . 
3 


= -V o /2+{V o2 /4+A(D/L)/pA j } h 


(E.2) 


(3.) 


Pump head is a function of losses and 


dynamic 




head. Assuming the losses AH^ to be 


only 




a function of jet velocity: 






AH. 

L, 


= Ka V. /2g 






H 


= \ - V o 2 /2g 


(E.3) 




therefore 






H 


= f(Vj) 






Q 


= f(Vj) 






giving 








SKP 


= f(Vj) 






SHP 


PgV j A j (K V . 2 /2g+h -V. 2 /2q) (E.4) 

550 n p n g 3 * ° 


(4.) 


Fuel v/eight burned for a constant SHP 


is : 




w F 


= SFC SHP t 


(E.5) 



where t is time in hours . 
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Appendix F. Program Flow Diagram 

CHART T lTie - pc )f f-o ()*(.<; 



I 

I 

I ‘i ;i r c o 
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• 07 
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I I 

I 1 P»NT r #, 1 

t 

I 

I 01 

f cpth ->pv / 

/ IPfAO / 

r v ! t fcrmat / 
r ioo / 

/ m~' THE L 1ST / 



I 

I NOT " 04 . 



* t J ST = V D ( ( I . 

* r? 25 F. 8 AnqL . 

* VJVf'J. NOT, 

* 1 FNGN 



I 

» 

I 

lM> U r rr*'STANTS 
I 
I 

I 0 

I = 3 . IMP 

I 

I **A • V». 

t 

I He - 12.0 

» 

I * • 32.17 

I 

| R*Ow * |.oa 



I 

I 

l r c T *?V=R O0*f, » 4 T 1 o 
ASSU^-ZJ T g£ 12 

i 

I 

I 06 

I r'>SGIl» a I 

» DISP/12.0 I 



» 

1 

I 

rnc&G 4T C8UISE 

PASE 3 X. t/C=L 4 

I 

I 

I 07 

I 3111 * I 

r iopicui i 

l /l*Hr**V'3(ll I 

{ *C*JV3 - 1.0(1 1 

I I 

I AJET * I 

I Qtl>/1 VJVO»Vj( 1 I ( I 



f 1 - 5 * 40(21 * I 

1 l.i *T 0*4G( l ) I 




I V3?2t * VDIIJO.* I 
I I 

I 0121 - I 

I * .1*7- IV 3 (2 1/2. ♦ I 

1 S<3Srr«V" 121 *V0( 2( 1 

1 ft-.O ♦ I 

I 7*0 40(2 1 I 

I /CJ>-*:w*4JFTJ ( » I 



( 

I 

| NOTF 0* 



fsPGT*' 'Y' L ~'D t> • 

2 i • (. y 



i 

i 



02 . 0 °» — >» 

4 I 10 



VJ( l( " OIK/ 4JET 

DfclHI l ,Q| = 

0 . 055 *VJ ( 1 1 * V J ( I ( 
/( 2. *G( 

OEtH(»,tO( = 

0. 0 1* VJ ( I (»VJ (II 
/ I 2 . 0 «G( 



I 

I 

2 • l l 



END Of Oi 
* L0JP7 • 




YES 



12 



• • N~> 

* FNI) °F 0* »-♦ 

• L - ,p 7 * I 

- - t 

• « I 

♦ res 

I . 07 . 

I .... 

I 

I 

) 

i w 



f r Ti IP ‘I T r >YST£H 
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05 /0?/ 7 ? 



•IT 



CHART TITLE - SUUft JUT TNI PU'^O 



02.1? >• 

PUMP JET "FS1GN 
program 

PUMP J r T Ppn^RAM 

SOLVES Thc pump, gear 
, AND FUEL WEIGHTS 
F CP 

BCTW T HP 
MUIT-PAR ALI EL 
CFNTP1FUGAL PIMP AND 
THF AXIAL PUMP 
THE PIMP DESIGN w|TH 
The l EAST 1DTAL 
WEIGHT IS G I VEN AS 
THE 

OPTIMUM DESIGN F"F 
the given input 

parameters 



PUMP CHAPCTEO 1ST K 
OATA 



K * CRUISE *»ciNT , J 
* TAKE (j p E P° INT 
PUMP OESIGN PCINT.J, 

I S AT T AKE OFF RASED 
ON TH0M4 »S CP ITER ION 
OFF DFS 1GN PC INT , K, 

1 S ANY OTHER Pul NT 
AL^NG THE QR A G TO 



" • <• 

* hsvi j i ,tr , 

• 0.5 DO * 



I TRUE 

I 



I 



e 1 * the value ~>n 

TH^MA* S CAV ITAT l r ’N 
IP I TFP I ON LESS Than 
TH r OIT 

^F PUNT, THD-l, THt 
I NP UCt A PLUS QN C 
AXIAL STAGE - IS 
REQUIRED 

I 

1 
I 



I XNSINSTG) «HPPt II I 
I • •D.TS/SOF T 1001 J I I 



I It p T UP 

f t: f 



l 



'IS 111 => 
l?*.0.« 001 Jl 

/ I «.S»',fOH|| J ] 

•II. “ 

3’ A T • » / l«JO"l J, II 
1 1 » • 1 1 . / I . I 

SHM l 1 
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/|P1»WPM| J, n 
• 01SI l l/t.0 .1"? 

I 



SH| instgi 
.GT. 0.41 



|FALS‘i . * . 

. 01 . 





SPEEO CUR 


1 


• 






1 


2 )1 1 




j 










1 

« 


I Tmom| , 0.055 


1 1 




1 


1 




1 


23 






1 


01 I 


1 


1 




1 


15 | 


1 


GGRAT INS TGI = 


1 












1 








1 


P?aF I 4, 1 ENGN1 


1 




1 


K = 1 


1 1 


1 


1 


1 


ETAPIJ.NSTG) = 


I 1 


1 


/RPHI J.NSIGI 


1 




1 




1 1 


• 0 3 


1 


1 


1.0 - 


1 I 


1 




1 




1 


J - 2 


1 1 


• • 


1 


1 


1 I3.TE5/D1SINSTG 


1 I 1 


1 


SHP(2,NSIS» = 


1 




1 




1 1 


• • 


I 


1 


1 • *0. 1 6 5M 1 • - 


1 1 


1 


Ph jw »G» 3 I 2 1 


1 




1 


1 PR NT = 6 


1 1 


• • 


TPJE 1 


1 


.0 1511 


1 1 


1 


•3PPI2I 


1 




1 




1 1 


• thomiji .lt. * 




1 




1 1 


1 /IS50.*XNGT*E TAP 


1 




1 


XPHMp = 


1 1 


■ TH3MI • 


1 1 


1 


OX = 03 IK I /00( J) 


1 1 


1 


(2, NS TGI*. ?R) 


1 




1 


FllATlNPUMP) 


1 1 


• • 


I 1 
















• - 


1 


1 


* • 

• . 


1 1 




1 

1 


1 

1 




1 








1 


1 


1 p AL SE . 


6 . 1 




1 


15 1 




1 


24 






1 note 


32 1 


1 


01 . 1 


• ■ 






. - - 














1 


... 1 


1 


fTAP |« ,11 = 


1 1 


1 1 


1 


H 




* 


PEG IN 00 LOOP 




I 


7 3 1 1 


1 


ET API J, 1 )• 1 PCI?, 


1 1 


10 


1 GEAR 


H 




* 


2 1 = J . 2 




1 


1 


1 


ll-Jx —2 . 


1 I 


1. 


1 


H 










1 


1 


1 


pci?, ? i*gx • 


1 1 


10 


1 


H 






1 


1 


1 


1 


1 


PC 12,31 1 


1 1 


1 3 


1 


M 










single stag inoucer 




















1 


03 1 


OESIGN 


1 




1 


1 




1 






• - 






1 


1 




1 


1 




1 






1 


OOUI => 


1 1 


1 


1 




1 


1 




1 


25 




I 


01 1 I/XPUMP 


1 1 


1 


1 


)PF OESIGN P "|NT pjmp | 










1 




1 1 


PUMP CHAR AIT E* 1 ST IC S 1 


p»"< 


1 


1 


xwGINSTGl = wG 


1 




1 


H SVI 1 I • HA - 


1 1 


1 


1 




1 


1 










1 


HE - DELHI 1,91 
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1 


1 




1 






1 


V)( | l*VOII 1 


1 1 


200 1 


ID 1 




1 


17 I 




1 






1 


/I 2.0*G1 


1 1 


— 






J 




1 N1TE 


25 






- 


1 


1 


1 1 


1 


PCI 1,21/PCI 1. 31 


1 1 


• 


BEGIN 00 L 3 OP 


• 






1 


1 


I XNSINSTG) = 149.5 


1 1 


1 


•0.5 


1 1 


• 


3 1 - 1 * 2 


• 






1 


0* 1 


» 


1 1 


1 




1 1 










* - 


— 




1 SMin * 


1 1 


1 


CB = 


1 1 




1 






1 


H P P | | ) = 


1 1 


1 XNSI l I 


1 1 


1 


PC (1 ,1 I/PCU.31 


1 1 


35 


.01 >1 






1 


1 VJI 1 1 *VJ (11- 


I 1 


1 •Il.o/THQMI J|) 


1 1 


1 




1 1 




1 


27 




1 
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1 1 


1 ••0.75 


1 I 


1 


HX ■ 


1 I 


• 








1 


/I2.0-G1 ♦ 


1 1 


1 


I 1 


1 


HPPI * (/HPP( J J 


I 1 


1 


ETAPP * 


1 




1 


OELHI 1,10) ♦ 


1 1 


1 PM I 111 - 0. 1 1 


1 1 









1 


PTAfM I.NSTGl 


1 




| 


OELHI 1 ,91 ♦ HF 


1 1 









1 


1 


1 


•ETAPI7.NSTG1 


1 




| 




1 1 


1 


1 




1 


1 


• - - 








1 


THOMII 1 * 


1 1 


1 


1 




1 


1 9 1 




1 






1 


HSVI 1 I /HP PI 1 I 


1 1 


1 


11 1 










1 






* - 





* | 


* — 




1 


R X = - C A • QX - 
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3 1 


2 B 






1 


1 
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1 1 


1 


SORT! }X»0<*ir 4« 


1 1 












1 


1 


1 
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1 


CA - CB! ♦ 
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1 


SHP(1 , N s T "» 1 = 
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1 


I 


1 CW = 347. I ■ XPUMP 
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1 


HX/PCI1 t3i 1 
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1 


RHJ4 *"j* 01 I 1 


1 






2 *05 


| 


| 


1 1 


1 




1 1 


1 


•HPPI I ) 


1 






• * 


I 


1 XXLP =1.70 


1 1 


1 


PPmIK.NSTGI = 


1 1 


1 


21 530. • XNGT« 


1 






• • 


1 


• 




1 


RPMI J , NSTG 1 •R X 


1 1 


1 


ETAPP 1 


1 




NO 


« • 


1 


1 


1 


1 




1 1 








♦ - 




• ENO OP DO 




1 


1 


1 


PhPAT = 0X/3X 


1 1 




1 








• i r)op-» • 


1 


1 


1 


* “ 








1 








• • 


1 


1 GN POINT P’JMo 


1 




1 


1 




1 








* • 


1 


SPPPD 


1 




1 

1 


1 

| 




1 t 








• 

1 YE S 


1 

1 


1 

1 


1 


DFY PUMP WEIGHT 


I 




t 5.31 








1 


1 


| 


12 1 




1 


1 












1 


1 


« — 






1 


1 














1 


1 PPM 1 J , NS T G 1 = 


1 1 




• 











< * )l 1ST ,1 = 

r-t- >1,1 nstgj*-?. j 

“LPINST^I - 
US|N< TGI * xxt P 

ipuu nst . i - 

0 . 7 ■* S« j 1SI TGI 

— 2*11. " 

SPAT AT ) 



-NTaei'io wATcP 
J ; | GMT 



I XWW(NSTG) * I 

I O.S>23«4PUP(N.> 1G) I 

I • °lp i ns r '» i i 

I •Hj.. u MPU«P | 

I 

I 

l 21 

I XNGT =• FlUATINGTl | 
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05/02/7? 



I'M •*. 



CHART T1TIF 



Sl*H Of lit I M, P(j«« O 



)•> 



0*. . 23 >* 

I 



01 



• FND JF 01 
• L n 0P? 



•3 



• .... 

I^S . (. . 

I . 27 . 

I .... 

I 

I 

I 



I 02 



1 1 1 


" 


H 


13 1 


fjfl 


H 


1. 1 




H 


10 1 




H 


11 1 




H 




1 

1 

1 


03 



I XW^INSTG) = *F I 



I 

I 

I 

TOTAL weight to 
DISPLACEMENT pa T t <7 

I 

I 

I 0<* 



wfat i nst g I =■ 

I XW01 NST 0 > • 
XWWINSTGI ♦ 
XWGINSTGI » 
XWFINSTGI I /DISP 



I 

I 

I 05 



/ WRITE TO DEV / 
/ 1PPNT / 

/ VIA FORMAT / 

/ 105 / 

/ FPCM THE l 1ST / 



I 

I 

I NOT E 0* 



list = vein* 

OISP. RANGE, 

V J V0 » NPUMP, not, 
IENGN 



I 

I 

I 07 

/ W«ITE to dev / 
/ IPPNT / 

/ VIA FOP* AT / 

/ 103 / 

/ F PGM THE L 1ST / 



I 

I 

I NOTE 0« 



l 1ST =» Oil*. * 

mppi i* , h s vi n, • 

thumi j I • 

I 

I 

I DO 



/ WRITE to dev / 
/ I PR N T / 

/ VIA FORMAT / 

/ 103 / 

/ <-PC« THE list / 



I 

1 

I NDTF 10 



• LIST * 31 21 • • 

• HPP12I, HCVI2I • 



1 II 

t watt- n n,v / 

/ |PSNT / 

/ VIA Cj»M 4 r / 

f IDA / 

/ snii the LIST / 



I 

I 

I N 1 T E 



12 
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I 

t 
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I 
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1 

I 

I 15 
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I 

I 
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LIST * 00 I J ) • 
0PM1J.NSTG1, 
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I 
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I 
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I 

I 
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I 

I 

1 71 
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K. = CRUISE POINT , J = TAKE OFF POINT 

PUMP DESIGN POINT, J, IS AT TAKE OFF BASED ON THOMA'S CRITERION 
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SHI ( 1) = G*HP'M J ) /( PI*RPM( J» 1 )*DIS( 1) /60. 1**2 
IF (SHI (NSTGI.GT .1.41 ) GO TC 701 

ETAP( J,NSTC) = 1. 0-C (3.666/DlS(NSTG) )**0. 1 65* ( 1 . - ,915 ) ) 
QX=GQ (K)/QO(J ) 

ETAP(K,1)=FTAP(J,1)*(PC(2,1 )*QX**2+PC(2 ,2 )*QX+PC (2,3)) 
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CALL FUEL 
XWF (NSTG) «UF 

TOTAL WEIGHT TC DISPLACEMENT RATIO 

WR A T ( NSTG) = (XWD (NSTG) + XWW ( NSTG ) +XW G ( NSTG ) +XW C ( NSTG ) ) /DI SP 
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